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A  Computer  Graphics  Package 


ABSTRACT 

1.  A  general  purpose  graphics  package  is  described.  The  package  is  an  updated 
version  of  a  graphics  library  first  developed  in  1969  on  an  ICL-1900  series 
computer  (described  in  AUVE  Acc  35647)  and  in  use  continuously  since  that  time 
with  various  modifications  and  informal  documentation  improvements.  It  is 
currently  implemented  on  a  DEC-MicroVAX  computer  and  various  hardware  devices 
including  the  HP  series  of  plotters  (using  HPGL),  Regis  devices  (VT125,  VT340), 
DEC  Workstation,  LN03  laser  printer,  and  certain  Ramtek,  Sigma  and  Tektronix 
devices.  - 

2.  Although  commercial  graphics  packages  are  available  this  one  is  simple  to 
understand  and  use,  and  together  with  its  familiarity  to  existing  users 
continues  to  provide  useful  graphics  support. 

3.  Examples  are  given  in  order  to  clarify  usage  of  the  software  and  lists  of 
subroutines  are  also  provided. 


1  BACKGROUND 


1.1  Recent  advances  in  computer  technology  and  software  have  dramatically 
changed  the  computing  scene.  Before  WIMPs  computers  were  used  mainly  by  skilled 
operators  or  in  rather  dedicated  applications.  Not  so  long  ago  suitable 
software  for  a  particular  application  was  unlikely  to  be  available  and  computer 
applications  relied  on  special  purpose  software  being  developed. 

1.2  In  scientific  research  printed  output  soon  acquired  a  low  reputation. 
Computer  solutions  to  complex  problems  frequently  produced  vast  quantities  of 
printed  numbers  and  data  was  transcribed  to  graphical  form  by  hand.  The  graph 
package  described  here  was  first  developed  in  about  1969  (ref:  AUVE  Acc  35647) 
for  use  with  a  Computer  Instrumentation  Ltd  incremental  plotter.  Under  computer 
control  this  machine  could  do  nothing  more  than  move  the  pen  one  increment  in 
one  of  eight  directions  or  raise  or  lower  the  pen.  Because  it  was  such  a  basic 
machine,  software  was  developed  so  that  scientific  applications  could  produce 
good  quality  graphics.  Features  included  font  designs  with  rounded  characters 
(unlike  the  poor  quality  angular  characters  which  were  available),  text  scaling 
and  rotation,  curve  fitting  to  allow  cusps  to  be  drawn  automatically,  and  indeed 
most  of  the  basic  features  currently  available.  The  package  has  served  the 
needs  of  scientific  graph  plotting  for  many  years  and  is  still  in  use,  including 
descendant  versions  adapted  by  various  users  at  ARE. 

1.3  Computer  graphics  has  become  a  very  important  and  sophisticated  subject  but 
has  many  different  applications,  including  computer  art  and  simulation.  Even 
scientific  applications  may  involve  complicated  graphics  which  are  not  suitable 
for  general  purpose  use.  Sophisticated  graph  packages  are  now  commercially 
available  and  have  the  advantage  of  providing  a  commercially  maintained 
standard,  such  as  GKS.  However,  this  package  continues  to  be  attractive  because 
it  is  simple  to  use,  it  is  familiar  and  can  be  easily  modified  or  extended  to 
cater  for  different  applications  or  new  hardware. 

1.4  This  report  marks  a  new  generation  of  the  software  for  compatibility  with 
the  window  environment  of  a  VAX  Workstation.  Some  changes  have  been 
incorporated  but  these  have  minor  implications  for  existing  programs  and  mainly 
deal  with  the  Workstation  aspects  and  additional  routines. 

2  INTRODUCTION 

2.1  ^In  scientific  research  graphical  output  from  computer  programs  is  essential 
and  can  be  easily  obtained  using  the  subroutines  described.  This  graphics 
package  provides  comprehensive  facilities  for  two  and  three  dimensional  work  and 
is  largely  device  independent.  The  subroutines  can  be  considered  in  three 
groups  as  follows: 

ier)  i  The  main  routines  which  actually  create  a  picture, 

(hj  'The  supporting  routines  which  modify  the  actual  behaviour  of  the  main 
routines ; 

(e^  Device  dependent  routines  which  are  not  normally  used  directly  by  the 
application  but  are  available  if  required  ' 

(d)  Special  routines  which  are  not  intended  for  general  use 


2.2  For  simplicity,  parameter  lists  have  been  kept  to  a  minimum,  and 
consequently  the  action  of  the  main  drawing  subroutines  depends  on  'hidden' 
parameter  values.  These  hidden  parameters  are  given  default  values,  but  may  be 
adjusted  by  means  of  the  supporting  subroutines.  For  example  the  subroutine 
PLOT_TEXT( ' text' ,N)  plots  text  using  the  default  character  size,  but  this  may  be 
changed  by  the  supporting  subroutine  CHAR_SIZE(Width, Height). 
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3  COORDINATE  SYSTEMS 


3.1  Assume  a  world  cartesian  coordinate  system  (X,Y)  is  defined  on  some 
imaginary  drawing  surface  such  as  an  A4  sheet  of  paper  with  its  origin  at  the 
lower  left  hand  corner  and  the  X-axis  horizontal.  The  units  of  both  X  and  Y  are 
mm.  A  local  cartesian  coordinate  system  (U,V)  may  be  defined  with  its  origin  at 
X=Xshift,  Y«Yshift  (see  ORIGIN)  and  scaled  such  that  one  unit  of  U  is  Xscale 
mm.,  and  similarly  for  Y  (see  SCALE).  Thus: 

X  =  U*Xscale+Xshif t  mm. 

Y  =  V*Yscale+Yshif t  mm. 


3.2  Initially  these  two  systems  (world  and  local)  are  identical  and  mapped  to 
the  actual  display  device  assuming  the  drawing  surface  is  A4  landscape  in  size 
and  orientation.  With  an  HP7750  plotter  for  example  the  mapping  is  to  true 
scale,  but  on  a  VT340  some  additional  scaling  is  necessary  in  order  to  represent 
the  A4  size  on  the  screen.  For  simple  applications  it  is  useful  to  imagine  the 
graphical  output  is  to  be  produced  on  an  A4  sheet  of  paper  in  landscape  or 
portrait  orientation.  However,  the  routine  PLOT_WINDOV  enables  the  world 
coordinate  system  to  be  modified  and  PL0T_VIEW  modifies  how  the  world 
coordinates  are  mapped  to  the  physical  display  device,  thus  enabling  large 
graphics  surfaces  to  be  created.  Workstation  graphics  may  create  up  to  4 
virtual  displays  with  4  windows  in  each.  Although  mapping  a  graph  to  a  display 
device  is  defined  in  terms  of  a  window  to  a  virtual  display  and  a  viewport  on 
the  device,  the  full  windows  environment  is  not  supported  on  some  devices. 

3.3  Valid  physical  graphics  devices  are  recognised  by  the  following  mnemonics: 

GF  A  pseudo  device  which  saves  a  file  copy  of  the  graphics  output. 

RG  A  Ramtek  series  9000 

VG  A  VT340  (or  VT125  with  some  restrictions) 

WG  A  Workstation  with  UIS  support 

HP  A  Hewlett  Packard  plotter  recognising  HPGL 

Any  hardware  combinations  can  specified  by  means  of  these  mnemonics  and 
additionally  the  mnemonic  TT  can  be  used  to  represent  either  VG  or  WG  determined 
by  the  login  device.  By  default  the  software  automatically  selects  the 
appropriate  display  device  you  are  logged  in  to  (if  it  is  recognised  as  a  valid 
graphics  device)  and  a  hard  copy  file  is  produced  (See  PLOTINIT) .  A  keyboard 
and  cursor  is  associated  with  the  first  interactive  device  selected. 

3.4  Special  facilities  exist  for  3D  plotting  and  are  described  in  greater 
detail  later. 

4  DEFAULT  CONDITIONS 

4.1  Initialisation  (see  PL0T_INIT)  resets  certain  parameters  to  default  values 
as  follows: 


ORIGIN 

0. 

0. 

SCALE 

1. 

1. 

CHAR  SIZE 

2.4 

3.0 

CHAR  ANGLE 

0.0 

CHAR  SLOPE 

0.0 

LINE  TYPE 

0 

0.0 

SELECT  PEN 

1 

HOVE 

0. 

0. 

Various  other  parameters  are  reset  by  device  specific  initialisation  routines 
(see  DEVICE  SPECIFIC  SUBROUTINES). 
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5  SUBROUTINE  DESCRIPTIONS 


5.1  Since  the  various  graphics  devices  have  different  capabilities,  it  is  not 
always  possible  to  ensure  that  general  purpose  routines  behave  identically  on 
the  different  devices.  In  particular,  character  drawing  will  be  subject  to 
resolution  differences.  Consequently,  it  may  be  necessary  to  make  some 
compromise  over  the  style  of  graph  for  a  given  device.  Alternatively  different 
devices  may  be  turned  off  and  on  to  allow  different  text  styles  to  be  used. 
However,  the  following  general  purpose  routines  provide  a  high  level  of  device 
independence,  where  the  parameter  type  specifications  are  INTEGER*4  (if  starting 
with  the  letters  I,J,K,L,M,N)  or  REAL*4  or  as  specified  in  the  description: 

A4_B0X(I) 

Draw  corner  marks  to  define  an  A4  box.  If  I>=0  the  long  side  is 
horizontal,  if  I<0  the  long  side  is  vertical. 

AXI  S_LIN( XO ,  YO ,  XI ,  Yl ,  SD ,  D ,  P  ) 

Draw  a  linear  type  axis  from  (X0,Y01  to  (XI, Yl),  with  distance  D  between 
tick  marks  which  are  of  height  P(mm.).  The  starting  position  through  the 
sequence  is  given  by  SD  expressed  as  a  fraction  of  D. 

AXIS_L0G(X0, YO, XI , Yl , SD, D, P) 

Draw  a  log  type  axis  from  (X0,Y0)  to  (XI, Yl),  with  distance  D  between 
cycles  using  tick  marks  which  are  of  height  P(mm.).  The  starting  position 
through  the  sequence  is  given  by  SD  expressed  as  a  fraction  of  D. 

CHAR_ABS ( CC V , CSH ) 

Reposition  at  the  character  coordinates  CSW*Character  width,  CSH*Character 
height,  relative  to  the  current  origin. 

CHARANGLE(PHI) 

Define  the  direction  in  which  characters  are  plotted  where  the  angle  PHI  is 
in  degrees  relative  to  the  x-axis  and  positive  anticlockwise.  (Default: 
0.). 


CHAR  P0SN(CSV,CSH) 

~See  CHARABS 

CHAR_REL( CSV , CSH ) 

'Reposition  at  the  character  coordinates  CSW*Character  width,  CSH*Character 
height,  relative  to  the  current  pen  position. 

CHAR_SIZE(Vidth, Height) 

Define  character  size  (mm.).  (Default:  2. 4, 3.0). 

CHARSLOPE(PHI) 

Causes  characters  to  be  slanted  at  an  angle  PHI  degrees  relative  to  the 
vertical,  with  positive  to  the  right.  (Default:  0.0). 

CURSOR (U,V, I char) 

Hove  the  cursor  under  operator  control  until  terminated  by  a  keyboard 
character  or  mouse  button.  The  cursor  coordinates  are  returned  in  U,V  and 
the  terminating  character  or  button  number  in  Ichar. 

CURVE (F,U_Min, U_Max , Tol ) 

Draw  a  smooth  curve  of  the  function  F  between  U_Min  and  U_Max.  Tol  defines 
the  smoothness  (segment  length  in  mm.). 

DRAV(U.V) 

Reposition  at  coordinates  (U,V)  with  pen  down  (see  plot).  Default 
coordinates  are  in  mm. 
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FILL_POLY(Xa , Ya , Na , Xb , Yb , Nb ) 

Using  the  currently  selected  colour,  fill  the  area  between  the  two  polygons 
defined  by  the  Na  points  in  arrays  Xa,Ya  and  the  Nb  points  in  the  arrays 
Xb,Yb,  with  polygon  a  circumscribing  polygon  b. 

FIT(U, V,N,M, SI , S2) 

Fit  a  smooth  curve  to  the  N  points  given  in  arrays  U,V.  If: 

(a)  M=0  the  slopes  are  not  defined  at  the  end  points,  and  will  be 
computed 

(b)  1  the  initial  slope  is  defined  in  SI 

(c)  2  the  final  slope  is  defined  in  S2 

(d)  3  both  slopes  are  defined  in  SI  and  S2. 

If  (U(N),V(N))=(U(1),V(1))  a  closed  curve  is  drawn. 

GREY(U,V,U) 

Draw  to  coordinates  U,V  using  intensity  or  colour  V. 

GREY_SCALE ( V_low , V  Jiigh ) 

Set  lower  and  upper  bounds  of  intensity  levels  for  V. 

GREYSET 

Initialise  grey  levels. 

GRID_LAT_LONG ( XO , YO , XI , Y1 , UO , VO , U1 , VI , IC ) 

Draw  a  latitude  v  longitude  grid  where  (X0,Y0)  (XI, Yl)  are  the  world 
coordinates  of  the  lower  left  and  upper  right  corners  of  the  grid 
respectively  (mm.)  and  (UO,VO)  (U1,V1)  define  limits  in  local  coordinates 
(degrees).  If: 

(a)  IC<0  only  the  scale  and  origin  are  computed. 

(b)  =0  as  above,  but  the  axes  are  also  drawn 

(c)  >0  as  above  with  the  axes  numbered. 

GRID_LIN_LIN ( XO , YO , XI , Yl , UO , VO , U1 , VI , IC ) 

Draw  a  linear  v  linear  grid  where  (X0,Y0)  (XI, Yl)  define  the  limits  in  mm., 
(U0,V0)  (U1,V1)  define  the  limits  in  local  units  and  IC  is  used  as  in 

GRID_LAT_LONG. 

GRID_LIN_LOG ( XO , YO , XI , Yl , UO , VO , U1 , VI , IC ) 

Draw  a  log  v  linear  grid  where  (X0,Y0)  (XI, Yl)  define  the  limits  in  mm., 
(U0,V0)  (U1,V1)  define  the  limits  in  local  units  and  IC  is  used  as  in 

GRID_LIN_LIN 

GRID_LOG_LIN(X0, Y0,X1 , Yl , UO, V0,U1 , VI, IC) 

Draw  a  linear  v  log  grid  where  (X0,Y0)  (XI, Yl)  define  the  limits  in  mm., 
(U0,V0)  (U1,V1)  define  limits  in  local  units  and  IC  is  used  as  in 

GRID_LIN_LIN 

GRIDLOGLOG ( XO , YO , XI , Yl , UO , VO , U1 , VI , IC ) 

Draw  a  log  v  log  grid  where  (X0,Y0)  (XI, Yl)  define  the  limits  in  mm., 
(U0,V0)  (U1,V1)  define  limits  in  local  units  and  IC  is  used  as  in 

GRID_LIN_LIN 

LINE_TYPE(I , P) 

Define  line  type  I  and  sequence  length  P  in  mm.  (Default:  0,0.). 

H0VE(U,  V) 

Reposition  at  coordinates  (U,V)  with  pen  up  (see  PL0T_ABS). 
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ORIGIN(X_Shif t.YShift) 

Redefine  the  local  coordinate  origin.  (Default:  0.,0.). 

PENDOVN 

Select  pen  down. 

PENUP 

Select  pen  up. 

PL0T_3D(X,Y,Z,IC) 

Plots  X,Y,Z  in  three  dimensional  coordinates,  as  if  viewed  from  an  azimuth 
relative  to  the  X-axis  and  an  elevation  from  the  X-Y  plane,  where  IC  is 
used  as  follows: 

(a)  IC=0  Initialise  using  azimuth  angle  -30.0  and  elevation  45.0 

(b)  =1  X,Y  are  the  azimuth  and  elevation  in  degrees 

(c)  -2  X, Y, Z  are  translation  factors 

(d)  =3  X,Y,Z  are  scale  factors 

(e)  =4  Move  to  X,Y,Z 

(f)  =5  Draw  to  X,Y,Z  without  hidden  line 

(g)  =6  Draw  to  X,Y,Z  with  hidden  line 

(h)  -7  Similar  to  6  ,ut  does  not  draw  (find  envelope) 

(i)  =8  Reset  envelope 

PL0T3DSURF ACE_XY (  X ,  Y ,  Z ,  N ,  M ) 

Draw  the  surface  defined  by  X,Y,Z(X,Y)  where: 

X  =  X(i)  (i=l,2. . .M) 

Y  =  Y(j)  ( j  =  1 , 2 . . ,N) 

Z  =  Z(i,j) 

Before  using  this  routine  a  3D  coordinate  system  must  be  initialised  by  use 
of  the  routine  PL0T_3D. 

PLOT  ABS(X,Y,IC) 

Reposition  at  world  coordinates  (X,Y)  expressed  in  mm.,  vhere  IC  is  used  as 
follows: 


(a)  IC<0  Use  pen  up 

(b)  =0  No  change  in  pen  condition 

(c)  >0  Use  pen  down 

PL0T_CHAR_1(M,N) 

Plot  the  character  number  N,  using  the  soft  character  set  number  M. 
PL0T_CHAR_  2 ( NDAT ) 

Plot  the  character  which  is  defined  by  data  in  the  INTEGER*4  array  NDAT. 
(See  source  code). 

PLOTCHARsymbol 

Plot  the  character  defined  by  symbol  as  follows: 


(a)  A  square-root  sign  when  symbol=SQR00T 

(b)  A  Greek  upper  case  character  when  symbol=name_U  and  name  is  one  of 
ALPHA,  BETA,  GAMMA,  DELTA,  EPSILON,  ZETA,  ETA,  THETA,  IOTA,  KAPPA, 
LAMBDA,  MU,  NU,  XI,  0MRIC0N,  PI,  RH0,  SIGMA,  TAU,  UPSIL0N,  PHI,  XI, 
PSI,  OMEGA 

(c)  A  Greek  lower  case  character  when  symbol=name_L  and  name  is  as 
before 
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These  symbols  have  been  provided  in  this  form  for  convenience  since  they 
commonly  occur  in  notation. 

PLOT_CONTOUR( ZO , Z , M , N , IC , SO ) 

Plot  the  contours  of  Z(X,Y)=Z0  where: 

X=i 

Y=j 

Z=Z(i|j) 

for  1=1... M,  j=l...N.  Special  effect  are  obtained  by  setting: 

IC=100*K+L 

where  integer  L  (<100)  is  the  pen  number  to  be  used  for  annotation,  and  if 
L  is  zero  no  annotation  is  produced.  The  integer  K  defines  whether 
interpolated  prints  on  the  grid  (i,j)  are  joined  by  straight  lines  (K=0)  or 
joined  by  a  smooth  curve  according  (K  non-zero).  If  selected  annotation  is 
placed  to  begin  at  a  distance  SO  along  the  contour  from  its  start. 

PLOT_DEVICE_(Text ) 

Parameter  Text  is  of  type  CHARACTERS *)  and  is  used  to  give  instructions  in 
the  form  keyword=value[ , . . . ]  where  [,...j  means  additional  expressions  may 
be  given.  The  form  of  each  expression  is  as  follows: 


Keyword 

Value 

DEVICE 

dd 

Where  dd  is  a  valid  device  name  in  use 

STATUS 

ON 

Enable  device  output 

OFF 

Disable  device  output 

ID-VD 

n 

Select  a  virtual  display  n  (1  to  4) 

ID-VD 

n 

Select  a  window  display  n  (1  to  4) 

TITLE 

"title" 

The  window  title 

no  device 

is  specified 

then  further  action  applies  to  all 

devices.  An  ID-VD  number  must  refer  to  a  virtual  display  previously 
created  by  PL0T_INIT.  An  ID-VD  number  may  refer  to  an  existing  window,  in 
which  case  it  is  selected,  or  to  a  new  window,  in  which  case  PLOTVINDOV 
must  be  used  to  create  the  required  window  (possible  preceded  by 
PL0T_VIEV).  If  no  virtual  display  or  window  is  specified  then  further 
action  applies  to  the  currently  sctive  window  (see  Examples).  An  empty 
character  string  (Text='  ')  causes  all  selected  devices  o  switch  status. 

PLOT_DEVICE 

Equivalent  to  PLOT_DEVICE_( '  '),  switch  status  of  currently  selected 
devices. 

PL0T_ERASE_(X0 , YO , XI , Y1 ) 

Erase  the  rectangle  defined  by  the  lower  left  corner  at  X0,Y0  and  the  upper 
right  corner  at  XI, Yl. 


PLOT_ERASE 

Erase  the  complete  display  or  viewport 


PLOT_FIN_(Text ) 

Parameter  Text  is  of  type  CHARACTERS*)  and  is  used  to  give  instructions  in 
the  form  keyword=value[ , .  . .  ]  where  [,...|  means  additional  expressions  may 
be  given.  The  form  of  each  expression  is  as  follows: 


Keyword  Value 

DEVICE  dd  Where  dd  is  a  valid  device  name  in  use 

STATUS  DELETE  Delete  the  file  copy 


7 


If  no  device  is  specified  then  further  action  applies  to  all  selected 
devices.  Normally  the  file  copy  is  saved  and  other  devices  are  deleted, 
but  the  file  copy  may  also  be  deleted  by  specifying  STATUS=DELETE. 

PLOTFIN 

Equivalent  to  PLOT  FIN_('  '),  save  the  file  copy  (if  selected)  and  delete 
other  devices. 


PLOT_INIT_(Text ) 

Parameter  Text  is  of  type  CHARACTER^*)  and  is  used  to  give  instructions  in 
the  form  keyword=value[ , . . .  ]  where  [,...]  means  additional  expressions  may 
be  given.  The  form  of  each  expression  is  as  follows: 


Keyword 

Value 

DEVICE 

dd 

Where  dd  is  a  valid  device  name  to  be  used 
namely  GF,  VG,  WG,  RG,  HP,  or  TT 

WINDOW 

FULL 

Select  and  map  to  the  full  screen  size  in  mm 

A4-H 

Select  and  map  to  A4  landscape 

A4-V 

Select  and  map  to  A4  portrait 

NONE 

No  window  created 

ID-VD 

n 

Create  virtual  display  number  n  (n=l,2,3,4) 

TITLE 

"text" 

The  window  title 

Any  hardware  combinations  can  specified  by  repeating  the  DEVICE=dd 
expression.  The  mnemonic  TT  can  be  used  to  represent  either  VG  or  WG 
determined  by  the  login  device.  If  no  device  is  specified  then 
'DEVICE=TT,DEVICE=GF'  is  assumed.  Use  of  TT  allows  an  interactive  device 
to  be  selected  without  GF.  If  no  window  is  specified  a  new  window  is 
created  and  mapped  to  A4-H.  The  window  option  NONE  is  used  to  over-ride 
the  default  window  creation.  In  this  case  no  window  is  created  at  this 
stage  (see  PLOT_WINDOW) .  An  empty  character  string  (Text='  ')  is 
equivalent  to  'WIND0W=A4-H' .  Initialise  the  hardware  and  software  as 
follows: 


(a)  Set  coordinate  scale  factors  (Xscale,Yscale)=(l,l) 

(b)  Set  coordinate  origin  (Xshif t,Yshif t)=(0,0) 

(c)  Set  character  size  (Xchsiz,Ychsiz)=(2.4,3.0) 

(d)  Set  character  direction  (PHI)  =(0.) 

(e)  Set  character  slope  (PHI)  =(0.) 


PLOT_INIT 

Equivalent  to  PLOT_INIT_('  ') 


PLOT_LIMITS( Umin , Vmin , Umax , Vraax ) 

Define  plotting  limits  in  current  plotter  units,  where  Umin, Vmin  defines 
the  lower  left  corner  of  a  rectangle  and  Umax,Vmax  the  upper  right  corner. 


PLOTMARK(I) 

Plot  an  identification  marker  at  the  current  position,  where  I  determines 
the  marker  shape,  in  the  range  1  to  9.  The  size  of  the  marker  is 
determined  by  the  current  horizontal  character  size  or  as  specified  by 
PLOTMARKSIZE(Size). 

P  LOT_MARK_SI ZE(Size) 

Determines  the  size  in  mm  of  marker  produced  by  PLOT_MARK.  Size=0.  resets 
to  the  current  horizontal  character  size. 


PLOTMODE(I) 

Define  the  writing  mode  as  follows: 
(a)  1.0  Replace  writing 
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(b)  1  Overlay  writing 

(c)  2  Complement  writing 

(d)  3  Erase  writing 

(e)  4  Negative  writing  off 

(f)  5  Negative  writing  on 

PLOT_NUM(Val,NrM) 

Plot  Val  (NB  Val  is  real)  as  a  character  string,  using  the  hard  character 
set,  with  N  figures  before  the  decimal  point  and  M  after.  If  M=0  an 
integer  format  is  used.  Positive  numbers  are  preceded  by  a  space,  and 
negative  numbers  by  a  -  sign. 

PLOT_NUM_SOFT ( Val , N , M ) 

Plot  Val  as  a  character  string,  using  the  soft  character  set  (see 
PLOT_NUM) . 

PLOT_REL(X,Y,IC) 

Reposition  at  relative  coordinates  (X,Y)  expressed  in  mm.  See  PLOT_ABS 
PLOT_TEXT(Text,N) 

Plot  the  N  characters  in  Text,  using  the  hard  character  set.  Special 

effects  can  be  obtained  as  follows,  where  n  is  an  integer  value: 

(a)  "nS"  Causes  n  spaces  to  be  output  (the  quotes  are  part  of  the 

syntax) 

(b)  "nC"  Causes  n  newlines,  relative  to  the  start  point 

(c)  "P"  Causes  the  character  ”  to  be  plotted. 

After  completion  the  pen  remains  at  the  end  of  the  current  line. 
PLOT_TEXT_SOFT(Text ,N) 

Plot  the  N  characters  in  Text,  using  the  soft  character  set.  Special 

effects  can  be  obtained  as  follows,  where  n  is  an  integer  value: 

(a)  "nS"  Causes  n  spaces  to  be  output  (the  quotes  are  part  of  the 

syntax) 

(b)  "nC"  Causes  n  newlines,  relative  to  the  start  point 

(c)  "nB"  Causes  characters  to  be  selected  from  set  n 

(d)  "P"  Causes  the  character  "  to  be  plotted. 

After  completion  the  pen  remains  at  the  end  of  the  current  line. 
PLOT_VIEV(Xmin ,  Ymin ,  Xmax ,  Ymax) 

Define  or  modify  a  viewport  for  the  current  window  using  Xmin,Ymin  as  the 
lower  left  corner  of  a  rectangle  and  Xmax, Ymax  as  the  upper  right. 

PLOT_WINDOV  (  Xmin ,  Ymin ,  Xmax ,  Ymax) 

Create  or  modify  a  window  to  the  current  virtual  display  using  Xmin, Ymin  as 
the  lower  left  corner  of  a  rectangle  and  Xmax, Ymax  as  the  upper  right. 
(See  PL0T_ DEVICE). 

READ_KB(Buf fer ,N1 ,N2) 

Read  data  from  the  currently  selected  keyboard  into  the  CHARACTER  string 
Buffer,  where  N1  is  the  maximum  number  of  characters  to  be  read  (unless 
input  is  terminated  by  <CR>,  and  N2  is  the  actual  byte  count  on  completion. 

READ_KB_NOECHO (Buffer, N1,N2) 

Similar  to  READJCB  but  without  echo. 

READ_KB_CHECK ( N ) 

Check  the  keyboard  input  buffer  and  set  N  equal  to  the  number  of  characters 
contained  in  it. 
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SCALE (X  Scale, Y  Scale) 

Redefine  tKe  coordinate  scale  factors.  (Default: 

SCFAC_LIN(H0,H1,SD,D) 

Calculates  suitable  scale  factors  for  axis  drawing,  where  HO, HI  are  the 
given  values  of  the  start  and  finish  points  and  D  is  the  calculated  tick 
length  interval  and  SD  the  position  of  the  first  tick  position  (as  a 
fraction  of  D).  D  and  SD  are  chosen  to  have  'nice'  values. 

SCFAC_LOG ( HO , HI , SD , D) 

Calculates  suitable  factors  for  log-axis  drawing,  where  D  is  the  cycle 
interval  and  SD  the  start  point. 


SELECTPEN(I) 

Select  pen  number  I.  (Default:  1). 


SELECTPENCOLOUR ( I pen , Red , Green , Blue) 

Define  the  shade/colour  of  pen  number  Ipen  as  specified  by  the  Red,  Green, 
Blue  values  (each  in  the  range  0.  to  1.)  with  special  effects  as  follows: 


Ipen 

0 

<0 


>=0 


Red  Green  Blue 

0.  0.  0. 

r  g  b 

r  g  b 


Reset  all  pens  to  a  standard  sequence 
Set  all  pens  to  shades  of  the  colour 
specified 

Set  pen  Ipen  to  the  colour  specified 


6  METHOD  OF  USE 

6.1  The  FORTRAN  examples  provided  demonstrate  the  use  of  many  of  the  routines 
described.  Each  program,  for  example  EXAMPLE_1,  should  be  compiled: 

FOR  Example_l 

linked  with  the  shareable  library  (named  GRAPH  here): 

LINK  Example_l , GRAPH/LIBR 
and  run : 

RUN  Example_l 

Example  1  uses  automatic  device  selection  and  therefore  a  file  copy  is  produced 
with  tKe  file  name  PLOT.PLT.  In  addition  a  display  copy  is  produced  if  the 
program  is  run  from  a  suitable  interactive  device.  The  file  copy  may  be 
displayed  on  a  suitable  interactive  device  by  use  of  the  command  DISP,  or 
plotted  by  the  command  PLOT.  Occasionally  it  is  also  necessary  to  link  with 
other  libraries  (such  as  MATHS).  Several  examples  are  listed  in  Annex  A. 
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Annex  A.  Examples 


A.l  The  following  examples  are  included  in  order  to  demonstrate  the  use  of 
various  routines.  Calculations  associated  with  developing  graphics  images  are 
not  proposed  as  realistic  or  optimum,  but  provide  a  simple  image  which  may  be 
displayed  or  plotted  using  the  routines  described. 
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3- JUL-3  9 


13:25 


EXAMPLE— 1 . FOR 


Paga  1-01 


c 

c 

c 

c 

c 

c 

c 

c 


c 

100 


PROGRAM  EXAMPLE_1 
Dascript ion : 

Graphics  sxaapls  showing  how  to  draw  •  siapla  graph.  Tha  dafault 
initialisation  is  usad  to  autoaiatically  dafina  an  A4  drawing 
surfaca  aappad  to  tha  davica  in  landscapa  oriantation. 


Author : 

C  Richardson,  ARE  (Portland) 


History: 


Xssua 


4  Novanbar  1983 


EXTERNAL  P1,P2 

Main  Entry  Point: 

PORMAT(A) 

CALL  PLOT_INIT 

CALL  CHAR_S IZE(1. 5,2.0) 

CALL  CHAR_SLOPE ( 20 . ) 

CALL  GRID_LIN_LIN( 60  . ,40.,260.,140.,0.,-1., 

CALL  MOVE ( 0 . ,0. ) 

CALL  DRAW( 360 .  ,0  .  ) 

CALL  CURVE { PI, 0  .  ,  360  .  ,1.  ) 

CALL  SELECT_PEN( 2 ) 

CALL  LINE_TYPE(2,10 .0) 

CALL  CURVE (T 2,0 .  ,  360.  ,1.  ) 

CALL  SELECTMEN  (  1  ) 

CALL  CHARTS I ZE( 1 . 8 , 2 . 4 ) 

CALL  PLOT_ABS (160. ,40. 1) 

CALL  CHAR_POSN<0 . ,-2 . ) 

CALL  PLOT_TEXT('X  ( Dag r aa s > '  , 1 1  ) 

CALL  PLOT_ABS (60. ,90. ,— 1) 

CALL  CHAR  POSR(-4.,0.) 

CALL  CHAR~ARGLE ( 90 . ) 

CALL  PLOT_TEXT ( ' S in ( X )  ',8) 

CALL  SELECTMEN  (  2  ) 

CALL  PLOT_TEXT ( ’ Cos ( X ) ' , 6 ) 

CALL  CHAR_ANGLE ( 0 . ) 

CALL  SELECT_PEN( 1 ) 

CALL  PLOT  ABS<20 . ,170 . ,-l> 

CALL  CHAR~AHGLE (-90.) 

CALL  PLOT_TEXT( 'Pig  l.',7) 

CALL  PLOTJTEXTC  A  Graphics  Exaapla',19) 
CALL  PLOT_TEXT('  Showing  Sin(x)  And  Cos(x)' 
ACCEPT  100, R 
CALL  PLOT  PIN 
STOP  'Finlshad' 

ERD 

PURCTIOR  PI ( X ) 

Fl-SIR(X/57. 29578 ) 

RETURN 

ERD 

PURCTIOR  F2(X) 

P2«COS(X/57. 29578) 

RETURN 

END 


IPunctions  to  plot 


llnitialisa  plottar 
iSalact  charactar  sisa 

360. ,1. ,1) 

(Draw  grid 

! Draw  a  lina  at  Y-0 . 

(Draw  sin(x) 

ISalact  pan  2 

ISalact  anothar  lina  typa 

I  Draw  cos ( x  ) 

ISalact  pan  1 

IPosition  for  X  labal 

IX  axis  labal 
IPosition  for  Y  labal 


IY  axis  labal 


IPosition  for  titla 


,26) 


lEnd  of  plotting 


■in ( x ) 


l  F2  is  cos ( x ) 
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Fig  Is  A  Graphics  Example  Showing  Sin  (x)  And  Cos  (x) 
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EXAMPLE-2  .FOR 


Paqa 


PROGRAM  EXAMPLE_2 

C  D«acr lpt Ion : 

C  An  axaapla  showing  a  function  which 

C  tha  routina  CURVE  and  nappad  to  tha 

C  oriantation. 

C  Autho  r : 

C  C  Richardson,  ARE  (Portland) 

c  History: 

C  Issua  1.0  4  Novanbar  1983 

EXTERNAL  Fn 

C  Main  Entry  Point: 

100  FORMAT (A) 


U0--10 . 

Ul-10 . 

V0  —  0.1 
Vl-1 . 1 

CALL  PLOT_INIT  ( ' WINDOW-A4-V » ) 

CALL  CHAR  SXZSTl .5,2.0) 

CALL  CHARTS  LOPE ( 2  0 .  ) 

CALL  GRID_LIN_LIN<  40.  ,50.  ,140.  ,170 .  ,U0,V0,U1 

CALL  CURVE ( Fn, UO, Ul, 1. ) 

CALL  CHAR_SIZE ( 1 . 8 , 2 . 4 ) 

CALL  MOVE ( (U0+Ul)/2. ,V0) 

CALL  CHAR_POSN ( 0 . , - 2 . ) 

CALL  PLOT  TEXT ( ' X ' , 1 ) 

CALL  MOVETU0 , (V0+V1 )/2 . ) 

CALL  CHAR_POSN ( — 4 . , 0 . ) 

CALL  PLOT_TEXT( 'y ' , 1 ) 

CALL  PLOT_ABS (  20  .  , 20  .  ,-l ) 

CALL  PLOT_TEXT( 'Fig  2. ',7) 

CALL  PLOT_TEXT('  A  Graphic*  Exampla',19) 

CALL  PLOT_TEXT ( '  Showing  a ',10) 

CALL  CHAR_POSN(0. ,0.25) 

CALL  CHAR  SIZE (1.5, 2.0) 

CALL  PLOT  TEXT ( ' - | X | ' , 4 ) 

ACCEPT  100, N 
CALL  PLOT  FIN 
STOP  'FinTshad' 

END 

FUNCTION  Fn ( X ) 

Fn»*xp(-ABS(X) ) 

RETURN 

END 


contains  a  cusp  plottad  by 
davica  in  portrait 


(Function  to  plot 


Ilnitialisa  plottar 
ISalact  charactar  siza 

,V1  ,1  ) 

(Draw  grid 
I  Draw  sin ( x  ) 

(Position  for  X  labal 

IX  axis  labal 
(Position  for  Y  labal 

JY  axis  labal 
(Position  for  titla 
(Plot  titla 


(End  of  plotting 


(Fn  is  axp(-x*x> 
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X 


Fig  3,  A  Graphics  Example  Shoeing  e  1x1 
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FILL-POLY-TEST . TOR 


Page  1-01 


C 

c 

c 

c 

c 

c 

c 

c 

c 

c 


c 

100 


PROGRAM  F1LL__P0LY_TEST 

P«8Ct iption : 

An  example  to  draw  closed  polygons  which  are  shaded  with  a  givan 
colour.  Tha  lattars  E  and  P  ara  drawn,  whara  tha  first  raquiras 
a  singla  polygon  to  ba  spacifiad  and  tha  sacond  raquiras  two. 
PLOT__INIT  spacifias  tha  portrait  format  A4-V. 

Author : 

C  Richardson,  ARE  (Portland) 


History  : 

Issua  1.0  6  Novanbar  1986 

Local  Variablas: 


REAL  X(60),Y(60),  ICoordinatas  for  lettar  P 

X 1 ( 1 3 ) /  ICoordinatas  for  lattar  E 

0.,  0  .  ,  1  5  .  ,  1 5  .  ,  5.,  5. ,15. ,25.,  5.,  5 . , 1 5 . , 1 5 . , 0 . /, 

Y1 ( 13  )/ 

0.  ,25.  ,25.  ,20.  ,20. ,15. ,15.  ,10.  ,10.  ,  5.,  5.,  0.,0./ 

Main  Entry  Point : 

FORMAT { A ) 


X(1)«0. 

Y(l)-0. 

X(2)-0 . 

Y ( 2 ) *25 . 

DO  n-0,24 
T-7 . 5*n 

X ( n+3 ) *7 . 5+7 . 5*SIND(T) 

Y(n  +  3)**17.5  +  7. 5*COSD(T ) 

ENDDO 

X(28)-5. 

Y(28)-10. 

X<  29  )«5. 

Y( 29)-0  . 

X(30>-5. 

Y( 30  > -20  . 

DO  n-0,24 
T-7 . 5*n 

X ( n+  3 1 ) -7 . 5+2 . 5*SIND(TJ 
Y(n+31 )-17 .5+2 . 5  *C0SD ( T ) 

ENDDO 
X ( 56 ) - 5  . 

Y ( 56 ) -1 5 . 

CALL  PLOT_INIT_< ' WINDOW-A4- V ' ) 
CALL  CHAR_SIZE( 1 . 5 , 2.0) 

CALL  CHAR  SLOPE (20.) 

CALL  SCALE ( 3 . , 3 . ) 

CALL  ORIGIN ( 50. , 50 . ) 

CALL  SELECT  PSNI4) 

CALL  F I LL_PO LY (XI, Yl, 13,0. ,0. ,0) 
CALL  SELECT_PEN( 2 ) 

CALL  MOVE ( XI (1)  ,Yl(l)  ) 

DO  n-2,13 

CALL  DRAW ( X 1 ( n ) , Y 1 ( n ) ) 

ENDD  ) 


IDefine  tha  outar  polygon 
!  for  lattar  P 


(Define  tha  inner  polygon 
I  for  lattar  P 


llnitialisa  with  portrait  format 


lEnlarga  tha  lattars  E,  P 

I  Position  of  lattar  E 
(Draw  lattar  E  shaded 
l  with  colour  4 
(Draw  tha  outline  using 
I  colour  2 


CALL  ORIGIN( 100 . , 50 . ) 

CALL  SELECTMEN  (  4  ) 

CALL  FILL  POLY(X( 1 )  , Y(  1  )  , 29  ,X(  30 )  , Y (  30  )  ,  27  ) 
CALL  SELECT_PEN( 2 ) 

CALL  MOVE ( X ( 1 )  , Y( 1  )  ) 

DO  n-2,29 

CALL  DRAW ( X ( n ) , Y (n ) ) 

ENDDO 

CALL  DRAW ( X ( 1 ) ,Y( 1 ) ) 

CALL  SELECT  PEN(3) 

CALL  MOVE(XT30 ) , Y( 30 ) ) 

DO  n-d,56 

CALL  DRAW ( X ( n )  ,Y(n)  ) 

ENDDO 

CALL  DRAW (X(30),Y(30)J 


(Position  of  lattar  P 
IDraw  lattar  P  shaded 
!  with  colour  4 
IDraw  the  outline  using 
l  colour  2 


IDraw  tha  inner  outline  using 
l  colour  3 


CALL  ORIGIN ( 0 .  , 0  .  ) 
CALL  SCALE(1 . ,1.) 
CALL  SELECT  PEN(l) 
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(Reset  scale  and  origin 
(Select  pen  1 
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FILL— POLX -TEST . FOR 


CALL  MOVE( 20 . , 20 . )  tPlot  a  title 

CALL  PLOT_TEXT( 'Fig  3. ',7) 

CALL  PLOT _ TEXT ( '  Example  Of  Drawing  A  Filled  Polygon',) 

ACCEPT  100, n 
CALL  PLOT_FIN 
END 
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Fig  3.  Example  Of  Drawing  A  FI  Hep  Polygon 
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PLOT-VAr-TEST . FOR 
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PROGRAM  PLOT_VAF_TEST 

C  Description : 

C  Example  using  variable  area  fill. 

C  Author: 

C  C  Richardson,  ARE  (Portland) 

C  History : 

C  Issue  1.0  6  November  1996 


C  Local  Variables : 

REAL  XX ( 201 ) , YY ( 201  )  ISimulated  data 

COMMON  T1,T2,T3,S,A,B,C,D 


C 

100 


Main  Entry  Point : 
FORMAT ( A ) 


N-201 
S-.  05 
A-100 
B»  2  0  . 
C-5  . 
D-2  . 


iNusber  of  points 
!Scale  factor 

{Data  generation  parameters 


CALL  PLOT  INXT  Unitialise  and  draw  a 

CALL  CHAR~SLOPE( 20 . )  l  rectangular  grid 

CALL  CHAR  SIZE ( 1 . 5 , 2 . 0  ) 

CALL  GRID_LIN_LIN{ 50.  ,50.  ,250.  ,186.  ,  0  .  ,  0  .  ,200.  ,17.  ,1) 

CALL  MOVeT1007, 0 . )  lAnnotate  x  axis 

CALL  CHAR_POSN(-5 . 5 , -2 . ) 

CALL  PLOT_TEXT( 'Time  (msec) ',11) 

CALL  MOVE ( 0 . , 8 . 5 )  lAnnotate  y  axis 

CALL  CHAR_POSN( -5 . , 0 . ) 

CALL  CHAR_ARGLE( 90 . ) 

CALL  CHAR_POSN(-8 . 5 , 0 .  ) 

CALL  PLOT_TEXT (’ Hydrophone  number ',17) 

CALL  CHAR_ANGLE( 0 . ) 

CALL  S CALE ( 1 . , 1 . )  IPlot  a  title 

CALL  ORIGIN ( 0 .  , 0  .  ) 

CALL  MOVE ( 20 . , 20 . ) 

CALL  PLOT_TEXT( ' Fig  4. ',7) 

CALL  PLOT_TEXT ( '  Example  Using  Variable  Area  Fill',) 


DO  >1,N  iCalculate  x  values 

XX( I )-I+49 . 

ENDDO 


DO  1-1,16 
Tl-20 . +5 .  *1 
T2-50 .♦7.5*1 
T3-80 .+10 . *1 
Y-8 . *1  +  50  . 

DO  IX— 1 , N 
X-IX 

YY ( IX ) -F ( X ) +Y 
ENDDO 

CALL  PLOT_VAF (XX , YY , N , Y+l  .  ) 
ENDDO 


t  Draw  trace  for  16  hydrophones 
ICalculate  times  for  3  events 
I  chosen  to  represent  some 
!  experimental  feature 

ICalculate  trace  for  this 
I  hydrophone 


(Draw  trace  with  filled  peaks 


ACCEPT  100,1 
CALL  PLOT_FXN 
END 
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FUNCTION  r<X) 
COMMON 


PLOT- VAF-TEST  FOR 


Tl  ,T2,T3,S,AyBrC,D 

XX-X-T1 

IF  (ABS(XX)  . LT .  0.000001)  THEN 

Y-A 

ELSE 

Y-A*SIN(XX)/(XX> 

ENDIF 

XX-X-T2 

IF  (ABS(XX)  .LT.  0.000001)  THEN 
Y-Y+A 
ELSE 

Y-Y+A*SIN(XX)/(XX) 

ENDIF 

XX-X-T3 

IF  (ABS(XX)  . LT .  0.000001)  THEN 

y-y+a 

ELSE 

Y-Y+A*SIN<XX)/(XX) 

ENDIF 

^s*rsiNix/5'>+c‘siN(x/2->"D#siN{x> 

RETURN 

END 
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PLOT- CONTOUR-TEST . TOR 


Paga  1-01 


PROGRAM  PLOT_CONTOUR_TEST 

C  Dascr iption : 

C  Exanpla  showing  a  contour  plot. 

C  Aut ho r : 

C  C  Richardson,  ARE  (Portland) 

C  History : 

C  Issua  1.0  6  November  1986 


C  Local  Variablas : 

REAL  2(15,11) 

C  Main  Entry  Point : 

100  FORMAT ( A ) 

101  FORMAT ( 16F8 . 3 ) 

M«15 
N-ll 
Ral3  . 

DO  i a 1  ,  M 
X-i-8  . 

DO  jal,N 
V-j-6 

Z0-R*  *2-X**2-Y*  *2 
IF  (20  .LE.  0.)  THEN 
Z(i, j)-0. 

ELSE 

Z ( i , j ) »SQRT ( Z 0  ) 

ENDIT 

ENDDO 

ENDDO 

CALL  PLOT_INIT  < ' WI NDOW=A4 -V ' ) 
CALL  CHAR_SIZE?1 . 5 , 2 . 0 ) 

CALL  CHAR_SLOPE( 20 . ) 

CALL  ORIGIN ( 10  .  , 50  .  ) 

CALL  SCALE ( 10 .  , 10  .  ) 

X*M 

Y*N 

CALL  MOVE ( 1  .  ,  1  .  1 
CALL  DRAW ( X , 1 . ) 

CALL  DRAW (  X  ,  Y  ) 

CALL  DRAW ( 1 . , Y ) 

CALL  DRAW ( 1 . , 1 . ) 

DO  J-2,N-1 
Y-J 

CALL  MOVE ( 1 . , Y) 

CALL  DRAW ( X , Y ) 

ENDDO 

Y«N 

DO  I a  2 , M— 1 
X*  I 

CALL  MOVE ( X , 1 .  ) 

CALL  DRAW ( X , Y ) 

ENDDO 
DO  i-l,M 
X-i 

CALL  MOVE (X , 1 . ) 

CALL  CHAR  ^OSW* -2 . 5 , -1 .  ) 

CALL  PLOT_NCM ( X- 8 .,2,0) 

ENDDO 
DO  j  - 1  ,  N 
Y-j 

CALL  MOVE ( 1 . , Y ) 

CALL  CHAR_POSN (-3.5,0. ) 

CALL  PLOT_NUM( Y-6 .,2,0) 

ENDDO 

CALL  MOVE (  (M+l .  )/2  .  ,  1  .  ) 

CALL  CHAR_POSN ( - 0  .  5  ,  - 2  .  ) 

CALL  PLOT_TEXT( 'X ' , 1 ) 

CALL  MOVE (l.,(N+l.)/r.) 

CALL  CHAR_POSN(-5  .  ,  0  .  ) 

CALL  PLOT_TEXT( ' Y '  , 1  ) 

CALL  SELBCT_PEN ( 2 ) 

IC-103 

SO-O  . 

DO  I2al2,19 
Z0aIZ/2  . 


(Grid  nap 


IMap  sisa 


tCalculata  Z  valuas 
I  for  tha  function 
!  X* *2+Y* *2+Z**2-R* *2 


! Ini t iali sa 

IChoosa  tha  scala  and  origin 

(Draw  outlina  of  grid 

1  Draw  grid 


lAnnotata  grid  X  axis 


lAnnotata  grid  Y  axis 


ISalact  anothar  colour 

(Draw  8  contours  with 
I  annotation 
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CALL  PLOT  CONTOUR ( ZO , Z ,M,N, IC , SO ) 

ENDDO 

CALL  SCALE ( 1 . ,1. ) 

CALL  ORIGIN { 0 . . 0 . ) 

CALL  MOVE (20. ,20. )  IPlot  a  titla 

CALL  SELECT  PEN(l) 

CALL  PLOT_T£XT{ 'Fig  5.',) 

CALL  PLOT_TEXT('  Contour  Plot  Of  Tha  runction  X',31) 

CALL  CHAR_POSN(0. ,0.25) 

CALL  PLOT_TEXT ( *  2 ' , 1 ) 

CALL  CHAR_POSN(0 . ,-0 . 25 ? 

CALL  PLOT^TEXT { *  +  Y',4) 

CALL  CHAR~POSN{ 0 .  ,0 . 25 ) 

CALL  PLOT _ TEXT ( ' 2 ' , 1 ) 

CALL  CHAR_POSN ( 0 . f -0 . 25 ) 

CALL  PLOT_TEXT('  +  Z ',4) 

CALL  CHAR_POSN ( 0  . ,0.25) 

CALL  PLOT_TEXT( ' 2 ' , 1 ) 

CALL  CHAR_POSN< 0 . ,-0 . 25 ) 

CALL  PLOT_TEXT<'  -  10 ',5) 

CALL  CHAR_POSN ( 0 . ,0.25) 

CALL  PLOT_TEXT( ' 2 ' , 1 ) 

CALL  CHAR_POSN( 0 . ,-0 . 25 ) 

ACCEPT  100,  IZ 
CALL  PLOT_riN 
END 
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25— SEP-89  16  ;  38 
PROGRAM 

C  Descript 

C 

C  Author : 

C 

C  History: 

C 


P LOT-SURF ACE-TEST.rOR 


PLOT_SURF ACE_TEST 

ion  : 

Example  showing  a  surface  plot. 

C  Richardson,  ARE  (Portland) 

Issue  1.0  5  December  1988 


Page  1-01 


C 


C 

100 


Local  Variables  : 


REAL  XX( 29 ) ,YY( 21 ) 

ZZ(29 ,21 ) 

Main  Entry  Point: 

FORMAT (A) 

MM- 2  9 
NN-21 
R-10  . 

DO  m-l,MM 

XX ( » ) -m- ( MM+1 . )/2 . 

DO  n-1 ,NN 

YY(n)=n-  NN+1 .  )/2  . 

ZO-R*  *  2 -XX  (  m  )  *  *  2-YY  (  n  )  *  *  2 
IF  ( Z0  .LE.  0 . )  THEN 
ZZ ( m , n ) -0 . 

ELSE 

ZZ ( a , n ) — SQRT ( Z0 ) 

ENDIF 

ENDDO 

ENDDO 


INode  coordinates 
{Surface  function 


{Surface  X  dimension 
{Surface  Y  dimension 
{Radius  of  sphere 
ICalculate  surface  height 
!  above  the  XY  plane 


{Allow  positive  values  only 


CALL  PLOT_INIT  {Initialise  graphics 

CALL  CHAR_S IZE(1. 5,2.0) 

CALL  CHAR_SLOPE( 20 . ) 

CALL  ORIGIN ( 50 .  ,  100  .  ) 

CALL  PLOT_3D( 0 . , 0 . , 0 . ,0 )  {Initialise  3D 

CALL  PLOT_3D( 5 . *MM/2 . . 5 . *NN/2 . , 0 . , 2 )  {3D  Origin 

CALL  PLOT  3D(5.  ,5.  ,5.  , 3 )  I  3D  Scale 


CALL  PLOT_3 D ( XX (  1),YY(  1),10.,4)  JDr 

CALL  PLOT_3 D ( XX (  l)rYY(  1),  0.r5) 

CALL  PLOT  3D ( XX ( MM ) , Y Y (  1),  0.,5) 

CALL  PLOT~3D(XX(MM) , YY(NN) ,  0.,5) 


CALL  CHAR_ANGLE (45.) 

CALL  CHAR__SLOPE  (45.) 

Y— 1 . 2*YY(I)-0 . 2*YY( 2 ) 

DO  m-1 , MM 

CALL  PLOT_3D(XX(m)  , YY ( 1  )  ,0.  ,4) 
CALL  PLOT_3D(XX(a) , Y, 0. ,5) 

CALL  PEN_UP 

CALL  CHAR_POSN (—4.0,0. ) 

IF  (MOD(m,2)  .EQ.  1)  CALL  PLOT 
ENDDO 

X-(XX( 1 )+XXfMM) )/2. 

CALL  PLOT_3D(X, Y , 0 . , 4 ) 

CALL  CHAR_POSN ( —3 .5,-1. ) 

CALL  CHAR_ANGLE ( -30 . ) 

CALL  CHARTS LOPE ( -30 . ) 

CALL  CHAR_POSN (-0.5,0. > 

CALL  PLOT_TEXT (  '  X  '  ,  1  ) 


{Annotate  X  axis 

tEnd  of  tick  mark 
{Plot  numbers 

(XX (a) ,2,0) 

lPlot  label  (X) 


X-0 . 8*XX(MM)+0 . 2*XX(MM-1 ) 

DO  n«l,NN  {Annotate  Y  axis 

CALL  PLOT_3D(XX(MM),YY(n) ,0. ,4) 

CALL  PLOT_3D(X,YY(n > ,0 . , 5 ) 

CALL  PEN_UP 

CALL  CHAR_POSN ( -0 .25,0.  ) 

IF  ( MOD ( n , 2 )  .EQ.  1)  CALL  PLOT  NUM ( YY ( n ) . 2 , 0 ) 

ENDDO 

Y-( YY( 1 )+YY(NN) )/2 . 

CALL  PLOT  3D(X,Y,0. ,4)  {Plot  label  (Y) 

CALL  CHAR~POSN( 3 . 5,-1 . ) 

CALL  CHAR_ANGLE( 45. ) 

CALL  CHAR~SLOPE (45.) 

CALL  CHAR~POSN(-0. 5, 0.  ) 
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CALL  PLOT_TEXT( •** , 1 ) 

CALL  CHAR__ANG  LE  (  -  3  0  .  ) 

CALL  CHAR_SLOPE( -30 . ) 

X»1.2*XX(1)-0.2*XX(2) 

DO  k«0,10 
Z-k 

CALL  PLOT  1D(XX(1) ,XX(1) ,Z,4) 

CALL  PLOT_3D(X,YY(i) ,Z,5I 
CALL  PEN  UP 

CALL  CHAR  POSH ( - 3 . 5 , 0 . ) 

IT  (MOD(X?2)  . EQ .  0)  CALL  PLOT_NUM ( Z , 2 , 0 ) 

ENDDO 

CALL  PLOT  3D(XX(1I ,YY(1 I ,5. ,4  I 
CALL  CHAR~*POSN  (-4.5,0.  ) 

CALL  PLOT~TEXT (  '  Z  '  ,  1  ) 

CALL  PLOT_3D_SURrACEJtY(XX,YY,ZZ,KM,NR)  IDr»w  surface 

CALL  CHAR_ANGLE ( 0 . )  lR.s.t  character  attributes 

CALL  CHAR  SLOPE (20.) 

CALL  ORIGINIO. .0. )  IR«*pt  origin  tnd  sc.lt 

CALL  SCALE ( 1 . , 1 . ) 

CALL  MOVE (20. ,20.  )  (Plot  titlo 

CALL  PLOT_TEXT( 'Pig  6. ',2) 

CALL  PLOT  TEXT ( '  surface  Plot  of  Tho  function  X',31> 

CALL  CHAR~POSN ( 0 . ,0.25) 

CALL  PLOT~TEXT  (  '  2  '  ,  1  ) 

CALL  CHAR_POSN (0 . ,-0 . 25  I 
CALL  PLOTJTEXTI'  +  Y’,4) 

CALL  CHAR_POSN(0. ,0.25) 

CALL  PLOT  TEXT I ' 2 ' , 1 ) 

CALL  CHAR“POSN( 0 . ,-0 . 25 ) 

CALL  PLOTJTEXTI'  +  Z',4) 

CALL  CHAR_POSN(0 . , 0 . 25) 

CALL  PLOT JTEXT (  ' 2 '  , 1 ) 

CALL  CHAR_POSN ( 0 . ,—0.25) 

CALL  PLOT_TEXT('  -  10', 5) 

CALL  CHARPOSNI 0 . ,0 . 25 ) 

CALL  PLOTJTEXTI ' 2 ' , 1 ) 

CALL  CHAR _ POSN (0.,— 0.25) 

CALL  SELECT_PEN( 2 ) 

ACCEPT  100, N 
CALL  PLOT_riN 
END 


{Annotate  Z  axis 

{End  of  tick  »ark 
[Plot  numbers 


1  Plot  label  IZ) 


> 
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CALL  PLOT  INIT  ( ' WIND0W-A4-V ' ) 
CALL  CHARTS 1ZE(3. 0,4.0) 

CALL  CHAR_SLOPE( 20 . ) 

CALL  MOVE (  20  .  ,150  .  ) 

CALL  PLOT_TEXT( 'Uppar  casa  ',11) 

CALL  PLOT_CHAR_ALPHA  U 

CALL  PLOT_CHAR_BETA_U 

CALL  PLOT_CHAR_GAMMA  U 

CALL  PLOT  CHAR_DELTA~U 

CALL  PLOT~CHAR_EPSILON_U 

CALL  PLOT_CHAR_ZETA  U 

CALL  PLOT  CHAR_ETA_U 

CALL  PLOT~CHAR_THETA  U 

CALL  PLOT_CHAR  IOTA_U 

CALL  PLOT  CHAR~KAPPA  U 

CALL  PLOT~CHAR  LAMDA~U 

CALL  PLOT_CHAR~MU_U  " 

CALL  PLOT_CHAR~NU_U 
CALL  PLOT_CHAR_XI_U 
CALL  PLOT_CHAR_OMRI CON  U 
CALL  PLOT_CHAR_PI__U  “ 

CALL  PLOT_CHAR_RHOJJ 
CALL  PLOT_CHAR_31GMA  U 
CALL  PLOT_CHAR_TAU_U~ 

CALL  PLOT  CHAR_UPSILON  U 
CALL  PLOT_CHAR_PHI  U 
CALL  PLOT_CHAR  XI  U 
CALL  PL0T_CHAR“PSI_U 
CALL  PLOT_CHAR_OMEGA  U 
CALL  MOVE (  20  .  ,  130  .  ) 

CALL  PLOT  TEXT ( ' Low* r  casa  ',11) 

CALL  PLOT~CHAR_ALPHA  L 

CALL  PLOT~CHAR_BETA_L 

CALL  PLOT_CHARJ3AMMA_L 

CALL  PLOT_CHAR~DELTA~L 

CALL  PLOT  CHAR  EPSILON_L 

CALL  PLOT”CHAR~ZETA  L 

CALL  PLOT~CHAR~ETA  L 

CALL  PLOT_CHAR~THETA_L 

CALL  PLOT  CHAR~IOTA_L 

CALL  PLOT~CHAR~KAPPA  L 

CALL  PL0T“CHAR“LAMDA“L 

CALL  PLOT_CHAR~MU_L  ~ 

CALL  PLOT  CHAR~NU~L 
CALL  PL0T“CHAR”XI  L 
CALL  PLOT  CHAR_OMRI CON  L 
CALL  PLOT~CHAR_PI  L 
CALL  PLOT_CHAR_RHO_L 
CALL  PLOT  CHAR_SXGMA  L 
CALL  PLOT~CHAR  TAU_L~ 

CALL  PLOT_CHAR”UPSTloN_L 
CALL  PLOT  CHAR-pHI_L 
CALL  PLOT~CHAR~XI  L 
CALL  PLOT~CHAR  PSl  L 
CALL  PLOT”cHAR~OHEGA  L 


CALL  CHAR  SI ZC ( 1 . 5 , 2 . 0 ) 

CALL  MOVETZO . , 20 . ) 

CALL  PLOT  TEXT (  ' Pig  7. ',7) 

CALL  PLOT”"TEXT (  '  Exaapla  Using  Oraak  cha r actara  '  ,  ) 
ACCEPT  107,1 
CALL  PLOT  PIN 
END 
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Upper  case  ABrAEZH® I KAMNX0TrP2TY$XPn 
Lower  case  apY6cS rje  i k Api>xonpdrv^x(pu 


Fig  7.  Example  Using  Greek  Charectera 
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PROGRAM  GREY_TEST 

C  Description : 

C  Example  plots  using  the  routine  GREY.  The  surfaces  which  are 

C  generated  do  not  represent  solutions  to  realistic  problems, 

C  but  provide  simple  examples. 

C  Author : 

C  C  Richardson,  ARE  (Portland) 

C  History : 

C  Issue  1.0  6  November  1986 


C  Local  Variables : 

INTEGER  ASKN 

C  Main  Entry  Point  : 

100  FORMAT (A) 

I=ASKN (' Pattern  number') 

CALL  PLOT  INIT 

CALL  CHARTS  XZE(1.8,2.4) 

CALL  ORIGIN ( 20 . , 0 . ) 

CALL  GREY_SET 

CALL  GREY  SCALE ( 0 10 . ) 

CALL  GREY_SELECT< 10 . ) 

CALL  MOVE ( 210 .  ,  160  .  ) 

CALL  CHAR_ANGLE (-90 . ) 

CALL  PLOT_TEXT( 'Example  Using  Grey 
CALL  PLOT_NUM( FLOAT < I ) , 2 , 0 ) 

IF  (I.NE.0)  GOTO  10 
DO  J=0,100 
G2-J/100 . 

Y-J*2 . 

CALL  MOVE ( 0 . ,Y) 

DO  1*0,100 
X=»  I  *  2  . 

Gl-I/100 . 

Z=5.  MG1  +  G2) 

CALL  GREY ( X , Y , Z ) 

CALL  READ_KB_CHECKf N) 

IF  (N.NE.0)  GOTO  91 
ENDDO 
ENDDO 
GOTO  90 

10  IF  (I.NE.l)  GOTO  20 

A* 1 . /2  0  0  0 . 

DO  J  —  50,50 
Y*  J  *  2  . 

CALL  MOVE ( 0 . , Y+105  ) 

DO  I*-50 , 50 
X»I *2  . 

Z-A* (X*X+Y*Y+X*Y) 

CALL  GREY ( X  +  100 . , Y+105 . , Z ) 

CALL  READ_KB_CHECK(N) 

IF  (N.NE.0)  GOTO  91 
ENDDO 
ENDDO 
GOTO  90 


(Five  different  patterns  are 
t  available  0-4 . 


{Initialise  grey  shades 
{Set  vertical  scale 
{Select  max  shade 
l  Plot  a  title 

Shading,  Pattern  Number', 42) 


{Pattern  0,  A  diagonal  ramp 
{Initial  grey  level  for  this  Y 

{Move  to  start 
{Scan  X  direction 

{Offset  grey  for  this  X 
{Actual  Grey  level  required 
{Draw  grey  segment 
{Test  keyboard,  Any  key  stops 
!  the  program 


{Pattern  1,  A  sqew  ellipsoid 


20  IT  (I.NE.2)  GOTO  30 

A-360./142.  {Pattern  2,  A  ripple 

C-.005 
DO  J*-50  ,  50 
Y* J  *  2  . 

CALL  MOVE ( 0 .  , Y+105.  ) 

DO  I--50,50 
X*I *  2 . 

R*SQRT (X*X+Y*Y) 

Z*10 .+5. *  COSD ( A  *R ) 

Z*Z  *  EXP ( -C  *  R ) 

CALL  GREY(X  +  100  Y  +  105.  ,Z) 

CALL  READ_KB_CHECK(N) 

IF  (N.NE.0)  GOTO  91 
ENDDO 
ENDDO 
GOTO  90 
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30 


40 


90 

91 


Ir  (X.NE.3)  GOTO  40 
A-COSD( 45 . ) *  3  60  ./X 4  2 . 

B«A 

C-0 .001 
DO  JaO.100 
Y-J*2. 

CALL  MOVE  f  0 ,  ,  y> 

DO  1-0,100 
X  — I  *  2  . 

R*SQRT(X*X+Y*Y1 

^*p?:^rA‘x+B‘i,+cosD,A*x-B**>> 

CALL  GREY (  X  ,  Y  ,  2 ) 

CALL  READ^KB  CHECK ( N ) 

If  (N.NE.i?)  GOTO  91 
ENDDO 
ENDDO 
GOTO  90 


°C  J-1,100 
Y-J*2  . 

CALL  MOVE ( 0 . ,Y) 

DC  1*1 , 100 

X— I  *  2 . 

T* AT AN 2 ( Y , X ) 

T«16 . • (T- . 785396  ) 
Z-l. 


*■10. •(.4+2| 

CALL  GREV(X,Y,2) 

CALL  READ_K8_CHECK(N) 
if  (n.ne.0)  goto  91 
ENDDO 
ENDDO 


ACCEPT  100,1 
CALL  PLOT  PIN 
END 


1  —  0  2 


JPafctarn  3,  a  vibrating  plata 


IP.tt.rn  4,  A  k,f  di.gr, » 
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Fig  8a.  Example  Using  Grey  Shading,  Pattern  0 


Fig  8b.  Example  Using  Grey  Shading,  Pattern  1 
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Pig  8c.  Example  Using  Grey  Shading,  Pattern  2 
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'*1*1 


Fig  8e.  Example  Using  Grey  Shading,  Pattern  4 
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PROGRAM  PLOT_SPHERE 

C  Description : 

C  Draw  a  sphere  on  a  chequered  board.  A  simple  calculation  is  used 

C  in  order  to  demonstrate  the  use  of  colour.  The  program  does  not 

C  model  perspective  drawing  or  light  reflection  in  a  sophisticated 

C  way . 

C  Author ; 

C  C  Richardson,  ARE  (Portland) 


C  History: 

C  Issue  1.0  26  January  1988 

C  Local  Variables: 


PARAMETER 

REAL 


CHARACTER  * 1 6 


Max«8  (Size  of  board 

X(0 :Max , 0 :Max ) ,  Y (0 :Max , 0 :Max ) , 

XP( 4 ) , YP( 4  )  , 

Xs(  0  :72 ,0:36)  ,Ys (0:72 ,0:36) 

Buf 


C 

C 

c 

c 


Main  Entry  Point : 

Calculate  various  coordinates  for  the  board  and  sphere.  The  board  is 
drawn  in  perspective,  but  the  sphere  Is  not  because  it  looks 
distorted  if  drawn  in  perspective. 

Max_Pen*l 50 
N-Max 
Xa-0  . 

Ya-220 . 

Xb-300  . 

Yb-220 . 

XO  *16  G  . 

Y0*40. 

Xl-140. 

Yl*190  . 

R*2 5  . 

CALL  INTERSECT l Xa ,Ya ,X0 ,Y0 , Xb ,  Yb,Xl , Y1 ,Xa0 , YaO ) 

CALL  INTERSECT ( Xa ,Ya , XI ,Y1 ,Xb, Yb ,X0 , Y0 , XbO ,Yb0 ) 

Sa*SQRT( (XO-XaO) **2+(Y0-Ya0 ) *  *2  )  (Some  constants  for 

Sa*Sa/ALOG ( 2 . )  I  drawing  the  board 

Sb-SQRT< (XO-XbO ) **2+(Y0— YbO 1**2) 

Sb«Sb/ALOG( 2 . ) 

SSa«SQRT( 1X0 -X a ) **2+(Y0-Ya ) *  *  2 ) 

SSb*SQRT(  (XO-Xb) **2+( YO-Yb ) *  *2 ) 

Aa*ATAN2D< Ya-YO,Xa-X0 ) 

Ab*ATAN2D( Yb-YO ,Xb-X0 ) 

X(0,0)*X0 
Y ( 0 ,0 )*Y0 
X(N,N)«X1 
Y ( N , N ) *Y1 

DO  1*0,  N 

U*TLOAT ( I )/N 
Db*Sb  * ALOG ( 1 . +U) 

XXb0»X0+Db*COSD( Ab ) 

YYb0%Y0+Db*SIND( Ab ) 

DO  J*0,N 

V-FLOAT( J)/N 
Da»Sa  * ALOG ( l .+V) 

XXaO*XQ+Oa*COSD< Aa ) 

YYa0*Y0+Da*SIND(Aa ) 

CALL  INTERSECT ( Xa ,Ya , XXbO , Y YbO , Xb , Yb.XXaO  ,  YYa0,X( I  ,  J  )  , Y(I , J)  ) 

ENDDO 

ENDDO 

DO  1*0,72  (Calculate  the  coordinates 

Th*5.*I  j  of  points  on  the  sphere 

DO  J*0 , 36 

Ph*5 .  *  J-  9 0 . 

U*R  *COSD ( Th ) *C0SD( Ph ) 

V*R*siND(Th ) *COSD( Ph ) 

W«R4R*SIND<5h) 

Xs ( I , J ) *X ( 4 , 4  )  *U 
Ys(I,J)*Y(4,4) *W 


(Corners  of  board 


(Calculate  all  the 
1  intersections  on 
!  the  board 


(Number  of  pens  to  use 
(Left  perspective  point 
(Right  perspective  point 
(Front  corner  of  board 
(Back  corner  of  board 
(Radius  of  sphere 
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ENDDO 

ENODO 

C  Using  th*  data  calculated  above  draw  the  board  and  sphere 

CALL  PL0T_INIT  (  '  DEVICE-TT  '  )  Unitialise 

CALI  CSLECT_PEN_CoLOUR(-l , 0 . ,1. , 0. )  ! Shades  of  green 

CALL  SELECTMEN  COLOUR  (  1 , 1  0  .  ,  0  .  )  !  Pen  1  set  to  Red 

CALL  SEL£CT_PEN~COLOUR( 2 , 0 . , 0 . , 1 . )  ! Pen  2  set  to  Blue 

CALL  SELECT_PEN  COLOUR ( Ma X_Pen , 1 .  , 1  1 .  )  [Last  Pen  set  to  white 

DO  Z«1,N  ~  J  Draw  the  board 

DO  J«1,N 

Xp ( 1 )«X( X-l , J-l )  [rind  the  four  corners 

YP( 1 )»Y< 1-1 ,  J-l J  !  of  a  square 

XP(2)*X(I,J-1) 

YP(2)-Y(X, J-l) 

XP ( 3 ) -X { X , J ) 

YP( 3 )*Y ( J , J ) 

XP{ 4 )»X< X-l , J) 

YP ( 4 )*Y ( J-l , J ) 

IP -MOD ( 1  +  J , 2 )  I  Alternate  the  colour 

CALL  SELECT  PEN(IP+1) 

CALL  riLL_POLY(XP, YP, 4 , , , 0 >  ’rill  the  square 

ENDDO 
ENDDO 

CALL  PLOT  MODE ( 0 ) 

Th0-300.  " 

PhO-45. 

Pi— COSD ( ThO ) *COSD( PhO ) 
p2»SIND|ThO) *COSD( PhO > 
p3-SIND(PhO ) 

DO  1-37,72 
Th-5 . *1 
DO  J-l  ,  36 
Ph«5. • J— 90  . 

XPU)»Xa(I-l,  J-l) 

YP(l)«Ys(X-l,  J-D 
XP( 2)-Xs(I, J-l ) 

YP ( 2  )-Ya ( X , J-l ) 

XP ( 3 ) -Xs ( I , J ) 

YP(3)«YsfI,J) 

XP ( 4 ) -Xs ( 1-1 , J ) 

YP(4 )-Ys (1-1 , J) 

Ql-COSO(Th) *  COS  D ( Ph 1 
<2  2-S  IND  (  Th  )  *  COSD  (  Ph  ) 

Q3— SIND ( Ph ) 

PQ-P1 *Q1+P2*g2+P3*Q3 
PQ-  (  PQ-vl  .  )/2  . 

IP-Max_P#n*PQ 

ir  (IP  .  GT .  Max_P*n)  IP-Max  Pen  .’Upper  limit 

IT  (IP  -LT.  3)  IP-3  ”  ! Lower  limit 

CALL  SELECT  PEN(IP) 

CALL  riLL_POLY(XP, YP, 4  ,  ,  ,  0 ) 

ENDDO 
ENDDO 

CALL  MOVE ( 50 . , 35 . ) 

CALL  CHAR  SIZE (3. 0,4.0) 

CALL  CHAR'S LOPE (20 . > 

CALL  $ELECT_PEN  (  Ma  x__Pen  ) 

CALL  PLOT_TEXT (’ Example  Showing  A  Shaded  Sphere',) 

CALL  REa0_KB_NOECHO ( Buf ,1,1) 

CALL  PLOT_MODE(l)  ’.Overlay  mode 

CALL  SELECT  PEN_COLOUR  (  0 , 0  .  ,  0  .  ,  0  .  >  .’Reset  colours 

CALL  PLOT _rlN 
END 


[Select  replace  writing 
[Set  the  lighting  direction 

trind  lighting  vector 
{Draw  sphere 

(Find  coordinates  of 
•  a  small  section 

[find  vector  normal 

•Find  angle  subtended 
[Convert  to  a  pen  number 
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PLOT-SPHERE. FOR 


SUBROUTINE  INTERSECT ( Xa 1 , Ya 1 , X* 2 ,  Ya2  , Xbl ,Ybl,Xb2,Yb2,X,Y) 

Sl-(Xal-Xa2)/(Yal-Ya2) 

S2* ( Xbl-Xb2 )/{ Ybl-Yb2 ) 

Y»(Sl*Ya2-S2*Yb2+Xb2-Xa2 )/(Sl-S2 ) 

X«Sl M Y-Ya2 )+Xa2 

RETURN 

END 


Fig  9.  Example  Showing  A  shaded  Sphere 
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PROGRAM  WINDOWS_TEST 

Description: 

An  example  to  demonstrate  the  usa  of  multiple  virtual  displays 
and/or  windows  on  a  workstation.  Two  virtual  displays  ara 
craatad  graan  and  blue  images  respectively.  In  each  virtual 
display  there  are  three  windows  each  mapped  to  the  screen  and 
selected  by  using  the  keys  A,B, 1,2,3  to  select  virtual  display 
1  or  2,  or  the  window  viewports  1,2,3. 

Author : 

C  Richardson,  ARE  (Portland) 


History : 

Issue  1.0  26  January  1988 

Local  Variables: 

REAL  Xs (0:72, 0:36), Ys(0:72, 0:36) 

CHARACTER  *  1 6  Buf 

CHARACTER*! 2  Ti tie ( 2 , 2 )/ 'Green  Sphe re ' , ' G r een  Cube 

'Blue  Sphere  ','Blue  Cube  '/ 


Main  Entry  Point: 

IP10»1 
IP11-100 
IP20-10 1 
IP21-200 

R-25  . 

DO  1*0,72 
Th-5. *1 
DO  J-0,36 
Ph»5 . * J-90 . 

U*R  *  COSD ( Th ) *COS  D ( Ph ) 
V-R*SIND(Th) *COSD ( Ph ) 
W*R+R*SIND(Ph) 

XS  < I , J ) *u 
Ys(I,J)»W 
ENDDO 
ENDDO 


(Pen  number  range  for  virtual 
!  display  number  1 
(Pen  number  range  for  virtual 
I  display  number  2 


(Calculate  the  coordinates 
i  of  points  on  the  sphere 


Initialise  virtual  display  number  1  and  create  window  number  1 
mapped  to  the  coordintes  0,0,  200,200 
CALL  PLOT  INIT  ( ' DEVI CE-WG , WI NDOW-NONE ’ ) 

CALL  PLOT~DLVICE  (' I D-WD» 1 , T ITLE* "Window  Number  1,1"') 

CALL  PLOT~VIEW(  0  .  ,0 .  ,  200  .  ,200  .  ) 

CALL  PLOT_WINDOW(  0  .,0.  ,200.  ,200.) 

CALL  CHAR  S I ZE ( 3 . 0 , 4 . 0 ) 

CALL  CHAR~SLOPE ( 20 . ) 

Draw  a  green  sphere  and  a  green  cube 
OO  I»IP10,IP11 

U-FLOAT( I-IP10 )/( IP11-IP10 ) 

CALL  SELECT_PEN  COLOUR ( I , 0 ., U , 0 . )  (Shades  of  green 

ENDDO 

CALL  SELECT_PEN  COLOUR ( IPH , 1 1 I . )  (Last  Pen  set  to  white 

CALL  ORIGIN( 50 .740 . ) 

CALL  SPHERE ( I P10 , IP1 1 , Xs , Ys )  (Draw  a  green  sphere 

CALL  MOVE (-25. ,-10. ) 

CALL  SELECT_PEN( IP11 ) 

CALL  PLOT  TEXT (Title(l,l),)  (Plot  title 

CALL  ORIGIN(150 . ,40 . ) 

CALL  CUBE (IP10,IP11)  (Draw  a  green  cube 

CALL  MOVE (-25.  ,-10.  ) 

CALL  SELECT_PEN( IP11 ) 

CALL  PLOT__T  EXT  (Title  (  2  , 1  )  ,  )  (Plot  title 

Initialise  virtual  display  number  2  and  create  window  number  1 

mapped  to  the  coordintes  100,50,  300,250 

CALL  PLOT_INIT_( ' DEVICE*WG , I D-VD* 2 , WINDOW-NONE ' ) 

CALL  PLOT_DEVICE  (' I D-WD-1 , TITLE* "Window  Number  2,1"') 

CALL  PLOT  V I E W ( 1  ?0  . ,50. ,300. ,250. ) 

CALL  PLOT“WINDOW(0 . ,0. ,200. ,200.) 

CALL  CHAR_S IZE(3. 0,4.0) 

CALL  CHARTS LOPE ( 20 . ) 

Draw  a  blue  sphere  and  a  blue  cube 
DO  I-IP20.IP21 

U*PLOAT( I-IP20 )/( IP21-I P20 ) 
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CALL  SELECT_PEN_COLOUR(I, 0.  ,0-  ,U)  l Shad* s  of  blu* 

ENDDO 

CALL  SELECT_PEN_COLOUR(IP21 , 1 . , 1 . , 1 . )  (Last  Pan  a*t  to  whit* 

CALL  ORIGIN (50 . ,40. ) 

CALL  SPHERE(IP20,IP21,Xs,Ys)  (Draw  •  blu*  spher* 

CALL  MOVE (—25. ,-10. ) 

CALL  SELECT  PEN ( IP2 1 ) 

CALL  PLOT_TEXT<Titl*(l,2> , )  (Plot  titl* 

CALL  ORIGINflSO. ,40. ) 

CALL  CUBEt IP20 , IP21 )  IDraw  a  blu*  cub* 

CALL  MOVE (-25. ,-10. J 
CALL  SELECT  PENUP21) 

CALL  PLOTJTEXT(Titl*(2,2) , )  ! Plot  titl* 

C  Using  virtual  display  1,  cr«*t*  window  2 

C  aapp*d  to  th*  coordintas  0,0,  100,100 

CALL  PLOT_DEVXCE  (' ID-VD-1 , ID-WD-2 , TITLE* "Window  Numb*r  1,2"*  I 
CALL  PLOT_VIEWf 0. ,0 . ,100. ,100. ) 

CALL  PLOT_WINDOW(0. ,0.  ,100  .  ,100.) 

C  Continuing  to  us*  virtual  display  1,  cr*at*  window  3 

C  mapped  to  th*  coordint*s  100,100,  200,200 

CALL  PbOT_DEVICE  I * ID-WD-3 , TITLE* "Window  Nu»b*r  1,3"*) 

CALL  PLOT_VIEW(  100. ,100.  ,200.  ,200.) 

CALL  PLOT_WINDOW(10  0  . ,0 .  ,  200  .  ,  100  ,  ) 

C  Using  virtual  display  2,  cr*at*  window  2 

C  mappad  to  th*  coordint*s  100,0,  200,100 

CALL  PLOT_DEVICE  (' ID-VO-2 , ID-WD-2 , TlTLE»"Window  Numbar  2,2"') 

CALL  PLOT  VIEWI100. ,50 . , 200 . ,150 . ) 

CALL  PLOT~WINDOW(0 . ,0 .  ,  100  .  ,  100  .  ) 

C  Continuing  to  us*  victual  display  2,  craat*  window  3 

C  mapped  to  th*  coordintas  200,150,  300,250 

CALL  PLOT_DEVICE_( *  ID-WD-3 , TITLE- "Window  Numbar  2,3"*) 

CALL  PLOT  VIEW(  200  .  , 150  .  ,  300  .  ,  250  .  ) 

CALL  PLOT”WINDOW(100. , 0 . ,  200 .  ,  100  .  ) 

C  Switch  batwaan  virtual  displays  1  and  2  by  pressing  ona  of  th*  keys 

C  A  or  B  and  switch  batwaan  windows  1,2,3  by  pressing  on*  of  th*  keys 

C  1,2,3.  Any  other  key  terminates  th*  program. 


Buf (1:1)-*1* 

DO  WHILE  <<<8uf(l:l)  .G£.  *  1  * ) . AND . ( Buf ( 1 : 1 )  . LT .  *4*)). OR. 

((Buf (1:1)  .GE.  'A'  )  .AND.  (Buf (1:1)  .  LT  .  *C'))) 

CALL  READ_KB  NOECHO ( Buf , 1 , I ) 

IP  (Buf (ill )“.EQ.  'A')  THEN 

CALL  PLOT_D£VICE  ('IO-VD-1')  (Salact  Virtual  display 

ELSE  IP  ( Buf (1:1)  TEQ.  'B')  THEN 

CALL  PLOT_DEVICE  ('ID-VO-2')  ISalaet  Virtual  display 

ELSE  IP  <  Buf  <  1  .*  1  )  .EQ.  '2')  THEN 

CALL  PLOT  DEVICE  ('ID-WD-2')  (Select  window  2 

ELSE  IT  ( Buf (1:1)  . EQ .  *3*)  THEN 

CALL  PLOT_DEVICE  ('ID-WD-3')  (Select  window  3 

ELSE 

CALL  PLOT_DEVICE  ('ID-WD-l')  (Select  window  1 

ENDIT 
ENDDO 


1 

2 


CALL  SELECT_PEN_COLOUR (0,0.,0.,0.) 

CALL  PLOT  TIN 

END 


(Reset  colours 
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Pafle  2-01 


SUBROUTINE  SPHERE ( IPO , IP1 , Xs , Ys  ) 

REAL  Xs  (0:72,0:  36) ,Ys( 0:72, 0:36)  , 

XP { 4 ) ,VPi4) 


Th0-300 . 

Ph0  =  4  5 . 

Pl-COSD(Th0)*COSD(PhO> 
P2=SIND(ThO )*COSD{PhO) 
P3-SIND(PhO) 

DO  1*37,72 
Th*5 . *1 
DO  J-1,36 

Ph«5 . * J-9  0 . 

XP<  1  )»Xs  ( 1-1 , J-l ) 
YP(l)*Ys(I— 1, J— 1 ) 

XP<  2 )*XS < I , J— 1 ) 

YP(2)»Ys (I, J-l ) 

XP( 3 )*Xs ( I, J) 

YP(3)*Y9(Ir J) 

XP ( 4 )*X9 { X-l , J ) 

YP(4  >«Ys(I-l, J) 

Ql-COSD(Th ) *  COS  D ( Ph ) 
Q2*SIND(Th) *C:SD( Ph ) 
Q3»SIND( Ph ) 
PQ-P1*Q1+P2*Q2+P3*Q3 
P<2»(PQ+1  .1/2. 

IP=PQMIP1-IP0  )+IP0 
ir  (IP  .GT,  IP1)  IP-IP1 
CALL  SELECT  PEN(IP) 

CALL  riLL_POLY(XP ,YP,4, , ,0) 
ENDDO 
ENDDO 
RETURN 
END 


! S • t  the  lighting  direction 
•Find  lighting  victor 

(Draw  sphere 

(Find  coordinetes  of 
!  a  small  section 


(Tind  vector  normal 

(Find  angle  subtended 

(Convert  to  a  pen  number 
f Uppe r  limit 


SUBROUTINE 

REAL 

R-25. 

XP(1)*0. 

YP(1 )»0 . 

XP ( 2 ) «XP ( 1 ) 
yp(2)*R 
XP(3)*-0,a66*R 
YP ( 3 ) ■! . 5  *R 
XP( 4 )=XP<3 ) 

YP( 4 )*0 . 5*R 

IP*NINT(  0 . 3333*  ( IP1-IP0  )  4-IP0  ) 

CALL  SBLECT_PEN( IP J 

CALL  FILL  POLY(XP,YP, 4, , ,0) 

XP( 3 )-0 . 8  66*R 

YP( 3 )-l  . 5*R 

XP(4)-XP(3) 

YP(4)*0.3*R 

IP*NINT( 0 .6667* ( ZP1-XP0 )+IP0 ) 
CALL  SELECT  PEN(IP) 

CALL  FILL  POLY(XP, YP, 4 , , , 0 ) 

XP(l)-0 

YP( 1 )*R 

XP(2>*0.866*R 

YP( 2 )-l . 5*R 

XP( 3  > -XP  < 1  ) 

YP ( 3 ) *2 . *R 
XP( 4 )*-0 .  $66*R 
YP(4)«YP(2) 

IP*IP1-1 

CALL  SELECT  PEN(IP) 

CALL  FILL  POLY ( XP , YP , 4 , , ,0) 

RETURN 

END 


CUBE( IPO , IP1  ) 
XP ( 4 ) , YP( 4 ) 


SUBROUTINE  INTERSECT ( Xa 1 ,Yal,Xa2,Ya2, Xbl ,Ybl ,Xb2 ,Yb2 ,X,Y) 

31*(Xal-Xa2  )/( Yal-Y«2 ) 

S2«(Xbl-Xb2)/(Ybl-Yb2) 

Y*(Sl*Y«2-S2*Yb2+Xb2-Xa2)/(Sl-S2) 

X*S1*(Y-Ya2)+Xa2 

RETURN 

END 
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Annex  B.  Summary  Of  Subroutines 


A4B0X(I) 

AXIS_LIN( XO , YO , XI , Y1 , SD , D , P ) 
AXIS_L0G(X0,Y0,X1,Y1,SD,D,P) 
CHAR_ABS(CSW,CSH) 
CHAR_ANGLE(PHI ) 
CHAK_POSN(CSV,CSH) 

CHAR  REL(CSV,CSH) 
CHAR-SIZE(Vidth, Height) 
CHARSLOPE(PHI) 

CURSOR ( U , V , Char ) 

CUR VE(F,Umin, Umax, TOL) 
DRAtf(U.V) 

FXLL_POLY( Xa , Ya , Na , Xb , Yb , Nb ) 

FIT(U,V,N,M,S1,S2) 

GREY(X,Y,Z) 

GREY  SCALE (Z  lov.Z  high) 

GREY  SET  “ 


GRID_LAT_LONG (X0,Y0,X1 ,ri,VO,VO,Ul,Vl,I) 
GRID_LIN_LIN(X0,Y0,X1,Y1,U0,V0,U1,V1,I) 
GRIDLIN_L0G<X0 , YO , XI , Y1 , UO , VO , U1 , VI , I ) 
GRIDL0G_LIN(X0 , YO , XI , Y1 , UO , VO , U1 , VI , I ) 
GRID_LOG_LOG(XO , YO , XI , Y1 , UO, VO, U1 , VI  I) 
LINE  TYPE(I,P)  ’  '  ' 

HOVETU.V) 

ORIGIN (X_Shift,Y  Shift) 

PEN  DOWN 
PENUP 


PLOT_3D(UU,W,W,IC) 

PLOT  3D  SURFACE(Z2, M, N, In terp) 

PLOT_ABS ( X , Y , IC ) 
PLOTCHARALPHA  L 
PLOTCHARALPHA  U 
PLOTCHARBETA  L 
P  LOTCHARBETAU 
PLOTCHARGAMMA  L 
PLOTCHARGAMMA  U 
P  LOTCHARDELTAL 
PLOTCHAR  DELTAU 
PLOT_CHAR~EPSILON  L 
PLOT_CHAR_EPSILON~U 
PLOTCHARZETA  L  ~ 
PLOTCHARZETA  U 
PLOT  CHARETA  L 
PLOTCHAR  ETAU 
PLOT_CHAR~THETA  L 
PLOT  CHAR~THETA-U 
PLOT  CHAR  IOTA  L 

plot“char~iota~u 

PLOT  CHAR  KAPPA  L 
PLOT~CHAR~KAPPA~U 
PLOT_CHAR_LAMDA~L 
PLOT_CHAP._LAMDA~U 
PLOTCHAR  MU  L  ~ 

PLOT_CHAR~MU~U 
PLOTCHAR  NUL 
PLOT  CHAR~NU~U 
PLOT_CHAR~XI~L 
PLOT_CBAR~XI~U 
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PLOT  CHAR_OMRICON  L 
PLOT'CHAR  OMRICON~U 
PLOT~CHAR"PI  L 

plot'char'pi  U 
plot"char"rho  L 

PLOT  CHAR'RHO  U 
PLOT^CHAR'SIGMA  L 
PLOT_CHAR  SIGMA_U 

plot_char'tau  L~ 
plot_char'tau"u 
plot_char"upsIlon  l 
plot_char~upsilon"u 
plot_char~phi  l  “ 
plotchar  phi-u 

PLOTCHAR  XI  L 
PLOT  CHAR_Xl~{] 

PLOT  CHARPSl  L 

PLOT~CHAR_PSI~U 

PLOT“CHAR  OMEGA  L 

PLOT_CHAR_OMEGA~U 

PLOTCHARSQROOT 

PLOT_CONTOUR( ZO , Z , M , N , IC , SO ) 

PLOT_COPY(FILE,N,M) 

PLOT  DEVICE 
PLOT'DEVICE  (Text) 

PLOT_ERASE  ~ 

PLOTJ2RASE  (XO, YO,Xl, Yl) 

PLOTFIN  ~ 

PLOT_FIN  (Text) 

PLOT~INIT 
PLOT JCNITJ  text) 

PLOTLIMITS (UmXn, Vmin, Umax , Vmax ) 
PLOTMARK(Select) 

PLOT  MARK  SIZE(Size) 

PLOT~MODEfr) 

PLOTJiUM(Val,N,M) 

PLOTJfUM  SOFT(Val,N,M) 
PLOT_RELTX,Y,IC) 

PLOT_TEXT ( Tex  t , N) 

PLOT  TEXT_SOFT(Text,N) 

PLOT  VAF(U, V,N, VO) 

PLOTVIEV(XMIN, YMIN, XMAX, YMAX) 

PLOT  WINDOW ( XMIN , YMIN , XMAX , YMAX ) 

READ_KB(Bu£fer ,N1 ,N2) 

READ  KB  CHECK(N) 

READ_KB~NOECHO (Buffer, N1 , N2) 
SCALE(X  Scale, Y  Scale) 

SCFAC  LIN(HO,Hl7DD,D) 

SCFAC  LOG(HO,Hl,DD,D) 

SELECf_PEN(I) 

SELECT_PEN_COLOUR (I, Red, Green, Blue) 


Annex  C. 


File  Contents 


A4B0X . FOR 

A4B0X(I) 


AXIS. FOR 

AXIS_LIN(X0,Y0,X1,Y1,SD,D,P) 
AXIS_LOG(XO , YO , XI , Y1 , SD, D, P ) 

CHAR_DATA. FOR 

CHAR_DATA ( Se  t , N , Char_bu  f ) 

COORDS. FOR 

SCALE (X  Scale, Y  Scale) 

ORIGIN  (  X__Shi  f  t ,  Y__Shi  f  t ) 

CURSOR. FOR 

CURSOR(X,Y,Char) 

CURVE. FOR 

CURVE(F,X_Min,X_Max,Tol) 


FIT. FOR 

FIT(X,Y,N,M,S1,S2) 

GRIDS. FOR 

SCFACLIN ( HO , HI , DD , D) 

SCFAC  LOG(HO,Hl , DD,D) 

GRID  LATLONG ( XO , YO , XI , Y 1 , UO , VO , U1 , VI , I ) 
GRID~LIN_LIN(XO , YO , XI , Y1 , UO , VO , U1 , VI , I ) 
GRID'LIN _LOG(XO , YO , XI , Y1 , UO , VO , U1 , VI , I ) 
GRID~LOG_LIN( XO , YO , XI , Y1 , UO , VO, U1 , VI , I ) 
GRID'LOG _LOG(XO , YO , XI , Y1 , UO , VO , U1 , VI , I ) 

PLOT  1 . FOR 

PLOT  INIT  (Text) 

'  PLOTINIT 

PLOTJfINDOtf  (Xmin ,  Ymin  ,Xmax,  Yraax) 

PLOT  VIEtf(Xmin, Ymin,Xniax) Ymax) 
PLOT~D£VICE_(Text ) 

PLOTDEVICE 
PLOT_FIN_(Text ) 

PLOT  FIN 
MOVE(U, V) 

DRAV(U.V) 

PLOTJREL ( X , Y , IC ) 

PLOT_ABS (X , Y, IC) 

PLOT2.FOR 

PEN  UP 
PEN~DOVN 
PLOT  MODE(I) 

PLOT~ERASE 

PLOT  ERASE  (XO,YO,Xl,Yl) 

SELECT  PEN~COLOUR( I , Red , Green , Blue ) 
SELECT~PEN?I) 

LINE_Ti?PE<I,P) 

CHAR  SIZE(VIDTH, HEIGHT) 

CHAR~SLOPE( PHI ) 

CHAR  ANGLE(PHI) 
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CHAR_SOFT 

CHAR_HARD 
CHAR_POSN ( CSV , CSH ) 

CHAR_REL( CSV , CSH ) 
CHAR_ABS( CSV , CSH ) 
PLOT_LIMITS ( Umi n , Vm in , Umax , Vmax ) 

PLOT  3D. FOR 

PL0T3D ( XX ,  YY , ZZ , CON ) 

PLOT  4. FOR 

FI LLPOLY ( Xa , Ya , Na , Xb , Yb , Nb ) 
GREY(U,V,tf) 

GREY  SET 

GREY~SCALE ( W_low , W_h igh ) 
PLOT_VAF(U,V,N,VO) 

PLOT  CONTOUR. FOR 

PLOTCONTOUR ( ZO , Z , M , N , IC , SO ) 

PLOTCHAR. FOR 

PLOT_CHAR_l(M,N) 

PLOT  CHAR_2(NDAT) 

PLOT_CHAR_2M ( NDAT , IC ) 
PLOT_CH(X,Y) 

GET_CH(NDAT , I , J ) 

SET  CH 


PLOTHARK. FOR 

PLOT_MARK( Select) 

PLOT_MARK_SIZE( Size ) 


PLOT_ SURFACE. FOR 

PLOT_3D_SURFACE ( ZZ , M , N , In t erp ) 

PLOTTEXTl . FOR 

PLOT_TEXT ( Tex  t , N ) 
PLOT_NUM(Val,N,M) 

PLOT_TEXT_2 . FOR 

PLOT_TEXT  SOFT(Text.N) 
PLOTNUMSOFT ( Val , N , M ) 

READJCB .  FOR 

READ_KB( Bu£ f er , N1 , N2 ) 
READ_KB_NOECHO( Buffer , N1 , N2 ) 
READ _KB_CHECK ( N ) 

SYMGREEK.FOR 

PLOTCHARALPHA 

PLOTCHAR  ALPHAL 
PLOT  CHAR~ALPHA_U 
PLOTCHAR  BETA 

PLOT_CHAR_BETA  L 
PLOT  CHAR_BETA~U 
PLOT  CHAR  GAMMA 

PEOT_CHAR  GAMMAL 
PLOT  CHAR-GAMMA_U 
PLOT_CHAR  DEETA 

PtOTCHAR  DELTAL 
PLOT  CHAR“DELTA  U 
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PLOT_CHAR_EPSILON 

PLOT_CHAR_EPSILON  L 
PLOT  CHAR  EPSILON  U 
PLOT_CHAR_ZETA 

PLOT_CHAR_ZETA  L 
PLOT  CHAR  ZETA  U 
PLOT  CHAR_ETA 

PLOT  CHAR_ETA  L 
PLOT'CHAR  ETA'U 
PLOT_CHAR_THETA  " 

PLOT  CHAR_THETA  L 
PLOT'CHAR  THETA'U 
PLOT  CHAR_IOTA  ~ 

PLOT_CHAR_IOTA  L 
PLOT  CHAR  IOTA  U 
PLOT  CHARKAPPA  ~ 

PLOT_CHAR_KAPPA  L 
PLOT  CHAR  KAPPA'U 
PLOT_CHAR_LAMDA  _ 

PLOT_CHAR_LAMDA  L 
PLOT  CHAR  LAUDA  U 
PLOTCHARMU 

PLOT  CHARMU  L 
PLOT-CHAR  MU  U 
PLOTCHARNU  ~  ~ 

PLOT  CHARNU  L 
PLOTCHAR  NU  U 
PLOT  CHARXI  _  “ 

PLOTCHARXIL 
PLOT  CHAR  XI  U 
PLOTCHAROMRICON-  ~ 

PLOTCHAROMRICON  L 
PLOT  CHAR  OMRICON  U 
PLOTCHARPI 

PLOTCHAR  PI  L 
PLOT  CHAR  PI  'U 
PLOTCHARRHO  ~  “ 
PLOTCHARRHO  L 
PLOT  CHAR  RHO'U 
PLOTCHARSIGMA  ~  ~ 

PLOTCHARSIGMA  L 
PLOT  CHAR  SIGMA~U 
PLOTCHARTAU 

PLOT  CHAR  TA'J  L 
PLOT  CHAR  TAU~U 
PLOT_CHAR_UPSILON" 

PLOT_CHAR  UPSILON  L 
PLOT  CHAR~UPSILON~U 
PLOTCHARPHI 

PLOT_CHAR  PHI  L 
PLOT  CHAR~PHI  U 
PLOTCHARXI 

PLOT_CHAR  XI  L 
PLOT  CHAR~XI  U 
PLOT_CHAR  PSI  ~  “ 
PLOT_CHAR  PSI  L 
PLOT  CHAR~PSI~U 
PLOT_CHAR  OMEGA  “ 

PEOT_CHAR  OMEGA  L 
PLOT_CHAR~OMEGA~U 


SYM  MATHS. FOR 

PLOT  CHAR  SQROOT 


